Decorrelation of the static and dynamic length scales in hard-sphere glass-formers 
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We show that in the equilibrium phase of glass-forming hard-sphere fluids in three dimensions, the 
static length scales tentatively associated with the dynamical slowdown and the dynamical length 
characterizing spatial heterogeneities in the dynamics unambiguously decorrelate. The former grow 
at a much slower rate than the latter when density increases. This observation is valid for the 
dynamical range that is accessible to computer simulations, which roughly corresponds to that of 
colloidal experiments. We also find that in this same range, no one-to-one correspondence between 
relaxation time and point-to-set correlation length exists. These results point to the coexistence of 
several relaxation mechanisms in the accessible dynamical regime of three-dimensional hard-sphere 
glass formers. 
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I. INTRODUCTION 

A recurring question about glass formation concerns 
the collective nature of the dynamics as one cools or 
compresses a liquid. If the phenomenon is collective, it 
should be characterized by the development of nontrivial 
correlations to which one or several typical length scales 
might be associated. One source of growing correlations 
has been clearly identified in connection with the increas- 
ingly heterogeneous character of the dynamics as the sys- 
tem becomes more sluggish. An associated length, com- 
monly referred to as "dynamical" , can then be extracted 
from multi-point space-time correlation functions [1-4]. 
In addition, several theories of the glass transition posit 
the existence of a growing "static" length accompany- 
ing a liquid's dynamical slowdown. This length is, how- 
ever, undetectable through standard measurements on 
pair density correlations, which have been shown to dis- 
play but unremarkable changes as the structural relax- 
ation slows down. 

Proposals for unveiling such a nontrivial static length 
include measures of the spatial extension of some locally 
preferred structure, as obtained from static correlations 
of a bond-orientational order parameter [5-15] or, via 
dimensional analysis, from the occurrence frequency of 
a given local arrangement [16-24]. Such proposals have 
long been advocated, but their usefulness remains uncer- 
tain [25, 26]. More recently, approaches that detect the 
growth in static correlations while staying clear of any 
specific proposal about local order, i.e., "order-agnostic" 
approaches, have been developed. Among these propos- 
als, we note patch repetition lengths [27, 28] , length scales 
extracted from information theoretic analysis [29, 30] or 
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from finite-size studies of the configurational entropy [31] , 
and other "point-to-set" correlation lengths [32-38] . 

Point-to-set correlations play a special role in the the- 
ory of the glass transition. They are more general than 
structural lengths based on a specific local-order descrip- 
tion, and are expected to provide upper bounds for the 
latter. As further discussed below, a type of point-to-set 
correlation length even enters in an upper bound for the 
relaxation time of the liquid [33] . Point-to-set correlation 
lengths can be studied by considering the distance over 
which boundary conditions imposed by pinning particles 
in a liquid configuration affect the equilibrium structure 
of the remaining (unpinned) particles. The original pro- 
posal, motivated by the random first-order transition the- 
ory [39], considered a cavity whose exterior is a frozen 
liquid configuration [32]. Other geometries of the set of 
pinned particles also allow one to extract point-to-set cor- 
relation lengths [35], although the lengths measured for 
different geometries need not coincide nor evolve in ex- 
actly the same way as temperature decreases or density 
increases [40]. 

It has also been suggested that a length scale could 
be obtained from the finite-size analysis of the relax- 
ation time itself, in a finite system with either periodic 
boundary conditions [31, 35] or with a pinned wall bound- 
ary [35, 38, 41]. Such a length has been called "dynam- 
ical" in Refs. 35 and 38, but it should be kept in mind 
that it is a priori different from the length characterizing 
the extent of the dynamical heterogeneity [42] . 

It is worth noting that, at present, none of the afore- 
mentioned lengths are directly accessible in experimental 
glass-forming liquids. The situation is slightly better for 
colloids and granular materials, but most of the informa- 
tion on these lengths must still be obtained from model, 
yet realistic, glass formers via computer simulations. 

Given this panorama of length scales that appear in 
the context of glass formation, many questions can be 
raised, among which the following two are central to the 
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glass problem. 

(i) Is the temperature (or pressure) evolution of these 
various lengths correlated? 

In other words, to which extent are the static lengths 
correlated amongst themselves, the dynamical lengths 
correlated amongst themselves, and the static and dy- 
namical lengths correlated with each other? Due to the 
limited growth of the static lengths that is generally ob- 
served, it seems hard to decide on the first issue. For 
the second, most data on dynamical lengths characterize 
the typical extent of the heterogeneous dynamics and are 
consistent with each other. As mentioned above, there is 
evidence that a "dynamical" length extracted from finite- 
size scaling behaves differently [38, 42], but we shall here 
mainly focus on the lengths extracted from four-point 
space-time correlation functions, which show no such de- 
viations. 

The third issue is more contentious. It has been force- 
fully advocated by Tanaka and coworkers [7-10] that the 
structural length extracted from the (static) correlations 
in a bond-orientational order parameter and the dynam- 
ical length obtained from four-point space-time correla- 
tion functions perfectly correlate. A somewhat differ- 
ent result was, however, obtained by Sausset et al. in 
their study of a glass-forming liquid on the hyperbolic 
plane [12], where the convergence of the two types of 
lengths is found to depend on the dynamical regime un- 
der consideration. A conclusion also at odds with that of 
Tanaka and coworkers has recently been reached by Xu 
et al. [15] for a polydisperse two-dimensional Lennard- 
Jones mixture and by Dunleavy et al. [30] for binary 
three-dimensional hard-sphere mixtures. 

(ii) Is the increase of the relaxation time due to the 
growth of any of the above lengths, or, with less com- 
pelling consequences, is it at least correlated to it? 

Guided by known forms of dynamical scaling near crit- 
ical points, several relations between relaxation time and 
lengths have been suggested and tested, such as a conven- 
tional power-law, r Q ~ £ 2 , and activated scaling behav- 
ior, log(r Q ) ~ with prefactors that possibly depend 
on temperature and pressure. Empirical correlations of 
one sort or the other have indeed been found with cither 
the conventional scaling and a dynamic length [43-47] 
(with z varying from 2 to 5) or the activated expres- 
sion and a dynamical [42, 48, 49] or a static length [7- 
9, 12, 23, 37] (with ip roughly varying from 1 to 2). 

Yet, as stressed by Harrowell and coworkers [50-52], 
correlation does not imply causation. The fundamental 
question to be addressed is therefore whether one can find 
a causal link between the increase of the relaxation time 
and that of any of the proposed length scales. An impor- 
tant result in this direction has been obtained by Monta- 
nari and Semerjian [33]: the relaxation time is bounded 
from above by an activated-like formula involving a static 
point-to-set correlation length £ps , 



r a < T exp (B £p S ) , 



(1) 



where tq sets the microscopic time scale and d 
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FIG. 1. Relaxation time (dotted line) and diffusivity (dashed 
line) of large spheres in two equimolar binary hard-sphere 
glass formers (7:5 squares; 6:5 circles) versus the reduced pres- 
sure pP/p, where p — N/V is the number density. Lines are 
guide for the eye. (Inset) Equation of state for the two mix- 
tures. 



the spatial dimension. The coefficient B depends on 
temperature and pressure, and is such that when £ps 
is about one particle size the right-hand side describes 
the "noncooperative dynamics" of the model [53]. 
According to the above equation, the relaxation time 
Tq, thus cannot diverge at a finite temperature nor at a 
finite pressure without the concomitant divergence of a 
static length. 

In this article, we address the above two questions by 
studying two glass-forming hard-sphere mixtures in d=3. 
We calculate a point-to-set length through the random 
pinning of a set of particles in an equilibrated configura- 
tion and a dynamical length via a four-point space-time 
correlation function. In conjunction with our previous in- 
vestigation of structural lengths associated with local or- 
der [25, 54], this study allows us to unambiguously show 
that the evolution with pressure (or packing fraction) of 
the static lengths decorrelate from that of the dynami- 
cal length in the range of relaxation times accessible to 
computer simulations. Additionally, no one-to-one corre- 
spondence between relaxation time and static length(s) 
is found. 



II. MODELS AND METHOD 

A. Model and simulation 

We consider two glass-forming hard-sphere systems 
in d=3; equimolar binary mixtures whose diameter ra- 
tios a\.02 are 7:5 and 6:5 (p\ sets the unit length) are 
selected to prevent crystallization. The properties of 
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these mixtures have been extensively characterized in 
Refs. [49, 55] and Refs. [56, 57], respectively. Equili- 
brated fluid configurations over a range of packing frac- 
tions ip with at least N = 1236 particles for the static 
length and N = 79, 104 for the dynamical length are ob- 
tained under periodic boundary conditions using a modi- 
fied version of the event-driven molecular dynamics code 
described in Refs. [26, 58]. Time is expressed in units of 
\J jima\ for particle of unit mass m at fixed unit inverse 
temperature /?. The diffusivity D is obtained by measur- 
ing the long-time behavior of the mean-squared displace- 
ment hm t ^ 00 ((Ar) 2 }(= -r«(0)] 2 ) = 2dDt, 
the pressure P is mechanically extracted from the colli- 
sion statistics, and the structural relaxation time r a is 
measured as explained in Sect. II C (Fig. 1). We now 
detail the procedures for extracting static and dynamical 
lengths. 



B. Static lengths 

We have already reported in Refs. [25 and 54] vari- 
ous measures of the spatial extent of the frustrated local 
tetrahedral order. To now compute a static point-to- 
set correlation length, we consider a system in which a 
fraction c of the particles of an equilibrium hard-sphere 
fluid configuration are pinned at random. Information 
on point-to-set correlations in the (bulk) fluid is then ob- 
tained from the long-time limit of the overlap between the 
original configuration and the configuration equilibrated 
in the presence of the pinned particles. If the reference 
and the final configurations are quite similar, the aver- 
age pinning spacing is shorter than the static correlation 
length, and the opposite if the two configurations are 
dissimilar. To measure the degree of similarity we use a 
microscopic overlap function [49] 



w mn {0,t) = e(a - \r n (t) - r m (0)|), 



(2) 



where a = 0.3a is chosen sufficiently small to enforce 
single occupancy for hard spheres. For a concentration c 
of pinned particles we therefore have 



Qc(t) = 



E 



y ' \m,n£B 



(0,t) 



(3) 



where the brackets denote an average over equilibrium 
configurations, the overline represents an average over 
the different ways to pin a fraction c of the particles of 
a given equilibrium configuration, and the sum is over 
all unpinned particles, with B denoting the set of pinned 
particles. 

By locating the rapid growth of Q c {oo) — Qo(°°) with 
c^ 1 / 3 , we extract the crossover length between small and 
large overlap. More specifically, we define the length £ p 
as the value of the average pinning distance for which 
the above overlap difference falls below 0.4. Note that 
the extracted length is not very sensitive to this choice, 



provided it is intermediate between low and high overlap. 
More methodological details can be found in Ref. [25]. 

In addition, we use Eq. (1) to estimate a lower bound 
on the growth of static point-to-set correlation lengths. 
From an Arrhenius-like argument for activation vol- 
umes [55], one expects B oc f3P for hard-sphere flu- 
ids [25, 54], indicating that the upper bound of r a di- 
verges with pressure even in the absence of any growing 
£ps, as when approaching T—0 for an Arrhcnius temper- 
ature dependence. In the low and moderate density fluid, 
the relaxation time indeed follows T a (P) ~ T a ,i ow (P) = 
To exp(K(3P) with K being a density-independent con- 
stant. One then finds that 



Cps(P) 
£ps,o 



> 



log[r a (P)/r ] 
log[r Q j ow (P)/r ] 



1/3 



(4) 



where £ps,o is the low-density limit of £ps and is related 
to K . The right-hand side of Eq. (4) thus provides a lower 
bound for the increase of any static length imposed by 
the dynamical slowdown. 



C. Dynamical length 

The dynamical relaxation of the fluid structure in the 
absence of pinning (c = 0) can also be obtained from the 
microscopic overlap function 




i(0,t) 



(5) 



which is similar to the self-intermediate scattering func- 
tion F s (q,t) for a wavevector q near the first peak of 
the structure factor S(q). The structural relaxation time 
T a can thus be approximated from the 1/e decay of 
F (t) [49]. 

Upon slowing down the fluid is known to exhibit fluc- 
tuations in particle mobility on a growing spatial range 
on the timescale r a [4] . The size of these regions defines 
a dynamical length £dyn, which can be extracted from 
the computation of a four-point space-time correlation 
function [49], 



G 4 (r;r Q 



V 



(N s (T a ))((N s (r a )) - 1) 



( ^2 w nn(0,T a )w mm (0,T a )S[r - r nm (0)] ) , (6) 
\n^m / 

where N s (T a ) — J2 n w nn{0, T a ) is the number of slow par- 
ticles on the structural relaxation timescale. The correla- 
tion length £dyn = ^4,(T a ) could be obtained from fitting 
G 4 (r;T a ) — Gi{r — > oo;r a ) to exp(— r/£ dyn )/r. This di- 
rect procedure is, however, numerically difficult because 
modulations arising from the fluid structure are superim- 
posed on the spatial decay of Gi{r,T a ). More robustly, 
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FIG. 2. For both the 6:5 (left) and the 7:5 (right) binary mixtures, the static (point-to-set) length £ p (dotted line) grows 
much slower than the dynamical length £dyn (dash-dotted line) relative to their value at ip, = 0.52, the onset of nontrivial 
dynamics [59]. The static length saturates and follows the bound given by Eq. (4) (dashed- line). The dynamical length £dyn is 
taken from Ref. 49 for the 7:5 mixture and is extracted from fitting the low-wavevector behavior of S4,(q,T a ) to Eq. (9) for the 
6:5 mixture. Lines are guide for the eye. (Inset) Collapse of the low-g behavior of S*4(g, T a ) using the Ornstein-Zernike form 
from Eq. (9) (solid line) for ip = 0.52 ; higher densities follow the form down to lower relative amplitudes. 



we use the Fourier space version of the function 
S 4 (<Z; r a ) = ^[(W(q;0,T a )W(-q;0,T a )) 



|(T4V(q;0,r a ))| 2 l, (7) 



where 



W(q;0,O 



w nrl (0,r ct )e 



-iq-r„(0) 



(8) 



which measures the structure factor of the regions in the 
system that remain immobile between times and r a . 
Fitting the low-g result to an Ornstein-Zernike form 



l + (g£dyn) 5 



(9) 



then provides the dynamical length £dyn- This analysis 
follows closely that of Ref. [49] for the 7:5 binary mixture, 
so we use these published values for this system, and only 
calculate the results for the 6:5 mixture. A good collapse 
of the low-q regime of 5*4 (q,T a ) is obtained by fitting 
Eq. (9) to the data (Fig. 2). 



III. RESULTS 

The results for the two glass-forming binary hard- 
sphere mixtures show that the point-to-set correlation 
length £ p increases only very modestly (by 80% for the 6:5 
mixture and less than 50% for the 7:5 mixture, see Fig. 2) 



for a density range over which the relaxation time r Q and 
the difiusivity D change by about 4 orders of magnitude 
(Fig. 1). As already shown in Ref. [25], the structural 
lengths associated with local order vary even less than 
£ p . Note that the bound £ps given by the right-hand 
side of Eq. (4) also increases only slowly and is compa- 
rable or less than £ p . Meanwhile, the dynamical length 
£dyn characterizing the spatial extent of heterogeneities in 
the dynamics markedly increases over the density range 
under study: £dyn grows by a factor of almost 4 for the 
6:5 mixture and 4.5 for the 7:5 mixture, with no sign of 
saturation (see also Ref. [49]). 

These results help answer the first question raised in 
the introduction. The growths of the dynamical and of 
the static lengths are not systematically correlated as the 
relaxation slows down. They even strongly decorrelate in 
the dynamical regime that is accessible to computer simu- 
lations, and to most colloidal experiments, which roughly 
corresponds to a 4 order-of-magnitude increase of the re- 
laxation time. The magnitude of this decoupling seems 
to be system dependent [60] . The divergence of the dy- 
namical and static lengths is more spectacular in the 7:5 
mixture than in the more weakly "frustrated" 6:5 mix- 
ture. For the latter, one could argue that the two types 
of lengths grow at the same pace at low pressure, while 
the relaxation time increases by, say, one order of magni- 
tude. Yet, even in this case, the two quantities eventually 
unambiguously part ways. 

We also investigate a possible correlation between re- 
laxation time and static length for the two hard-sphere 
binary mixtures. The outcome displayed in Fig. 3 shows 
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FIG. 3. The increase of the logarithm of the relaxation time 
log[r a (P)/ro] with f3P(, p for two binary hard-sphere mixtures 
(7:5 squares; 6:5 circles). The lines are guides for the eyes 
that merely illustrate the nonlinearity of the dependence. In- 
set: log[r a (P)/r*](/3P//3P*) versus £p/C* where time, pres- 
sure, and length are normalized by their value at y>* = 0.55. 
The lowest density point of the main figure is then dropped. 
Note that because of the implicit offset in the intercept, the 
roughly linear appearance of the plots does not have a simple 
interpretation in terms of an activated scaling expression with 



log[r Q (P)/ro] versus /3P£ P . The observed behavior is 
quite different from that reported in Ref. [37], where 
a data collapse for all three studied three-dimensional 
glass-forming liquids was obtained when using a sim- 
ple linear dependence corresponding to an activated-like 
scaling expression with ip = 1. Here, we find that a 
linear fit clearly does not describe the data and that 
the two hard-sphere mixtures cannot be collapsed onto 
a unique master curve. As seen in the inset of Fig. 3, 
this is true even when restricting the analysis to densi- 
ties above <p = 0.55. (Note that the "onset" value above 
which nontrivial glassy dynamics is reported to be around 
0.52 [59].) It should be stressed that the span of relax- 
ation times described in Ref. [37] is rather limited, cov- 
ering only 1.5 orders of magnitude. 

This second (empirical) finding of a nonuniversal re- 
lation between time and length, together with the very 
modest increase of all static lengths potentially associ- 
ated with the slowdown of relaxation, casts doubts on 
the existence of a general, one-to-one, causal relation be- 
tween the two quantities. In the simulation accessible 
regime of three-dimensional hard-sphere glass formers, 
at least, this relation is not observed. 



IV. CONCLUSION 

By studying two different three-dimensional binary 
hard-sphere glass formers, we have shown that the dy- 



namical length associated with the increasing heteroge- 
neous character of the dynamics and the various static 
lengths that have been put forward to explain the collec- 
tive nature of the slowdown of relaxation, be they struc- 
tural lengths associated with local order or point-to-set 
correlation ones, unambiguously decorrelate as density 
increases. This result is obtained in the dynamical regime 
that is accessible to computer simulations, which covers 
4 orders of magnitude in relaxation time and diffusivity. 
All considered lengths increase with density and relax- 
ation time, but the dynamical length grows much more 
rapidly than the static ones. 

This finding is sufficient to rule out a general princi- 
ple tying together the evolutions of dynamical and static 
lengths in glass-forming systems. It is possible however 
that the absence or presence of correlation between the 
quantities depends on the dynamical regime under con- 
sideration as well as on the type of material. A strong 
correlation among length scales is expected if glass for- 
mers are close enough to a putative thermodynamic crit- 
ical point, whether avoided or unreachable [61]. This 
phenomenon is what is predicted for instance in weakly 
frustrated systems, as is possibly observed in some two- 
dimensional systems [7-9, 11-13, 15], and near a random 
first-order transition to an ideal glass [62] . 

The link between relaxation time and static correla- 
tion length derived by Montanari and Scmerdjian [33], 
which we have somewhat heuristically extended and used 
in this paper, puts a bound on the contribution that can 
be attributed to a collective or "cooperative" activated 
mechanism driven by the growth of a static length scale. 
We have seen that this contribution stays rather modest 
in dynamical the range studied. In conjunction with the 
fact that no master curve is found to collapse the depen- 
dence of the relaxation time on the static length for the 
two different hard-sphere mixtures, this result points to 
the absence of a one-to-one correspondence between time 
and static length in the range under study, and therefore 
to the absence of a direct causal link. Lengths and time 
all grow in concert but it is impossible to assign them, 
on the basis of simulation data alone, a unique origin. A 
couple of factors indeed blur this issue: 

• Several relaxation mechanisms are likely to coexist 
and entangle in the accessible regime, with con- 
tributions coming from both cooperativity-driven 
effects characterized by a static length and facili- 
tation or flow/mode-coupling effects characterized 
by a dynamical length. 

• In a regime where length scales are modest (espe- 
cially the static ones) , various determination meth- 
ods may lead to different results, which further 
obscures the search for a causal relation with the 
structural relaxation slowdown. 

It is possible that there exists an asymptotic regime 
in which one mechanism dominates and characteristic 
lengths become very large so that one can causally at- 
tribute the slowdown of relaxation to the growth of a 
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unique typical length associated with the extent of col- 
lective/cooperative behavior driving the dynamics. This 
phenomenon is what several theories that involve a singu- 
larity of one form or another predict. Yet, such a regime, 
because of timescale limitation or irreducible frustration 
or obstacles, appears to be out of reach of (present day) 
simulation studies. Changing the curvature of space to 
increase the static length [11, 12] or increasing the di- 
mensionality of space to decrease it in the dynamically 
accessible regime [26, 54, 63-66], may therefore be more 



productive numerical approaches to understanding the 
glass problem. 
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